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We obtain the effective action of tunnel-coupied quantum dots, by modeiing the system as a Lut- 
tinger iiquid with muitipie barriers. For a doubie dot system, we find that the resonance conditions 
for perfect conductance form a hexagon in the plane of the two gate voltages controlling the density 
of electrons in each dot. We also explicitly obtain the functional dependence of the interaction en- 
ergy and peak-splitting on the gate voltage controlling tunneling between the dots and their charging 
energies. Our results are in good agreement with recent experimental results, from which we obtain 
the Luttinger interaction parameter K = 0.74. 



PACS number: 71.10.Pm, 73.23.Hk, 73.63.Kv 

Coulomb blockade CB [1] in quantum dots has gained 
in importance in recent years, due to its potential for 
application as single electron transistors. Recent exper- 
iments [2-4] on coupled quantum dots have revealed a 
variety of features, the most salient of which is the con- 
ductance peak splitting controlled by inter-dot tunneling. 
The peaks split into two for a double dot system and into 
three for a triple dot system. 

Earlier theoretical studies [5] on coupled dots studied 
the effect of tunneling and inter-dot capacitances on the 
CB, but only a few [6-8] focused on junctions with only 
one or two channels. Motivated by a) the double dot 
experiments, b) the analysis of Matveev [9], which maps 
a quantum dot formed by a narrow constriction that al- 
lows only one transverse channel to enter the dot, to a 
one-dimensional wire model and c) recent numerical ev- 
idence [10] that justifies modeling of quantum dots by 
Luttinger liquids, we model the coupled dot system as a 
one-dimensional Luttinger liquid (LL) with large barriers 
and study the effective action in the coherent limit. 

We obtain a simple expression for the interaction en- 
ergy in terms of the charging energies of the two dots and 
the gate voltage controlling the inter-dot tunneling, and 
find that tunneling between the dots splits the conduc- 
tance maxima. With (without) inter-dot tunneling, the 
conductance (resonance) maxima, for a double dot sys- 
tem, form a rectangular (hexagonal) grid in the plane of 
the two gate voltages controlling the density of the elec- 
trons. The peak splitting grows to a maximum when the 
two dots merge into a single dot. From the experiments, 
we estimate the interaction energy and also the values 
for the Luttinger parameters. 

Although a classical analysis with an inter-dot capac- 
itance reproduces the experimental results qualitatively, 
the value of the required inter-dot capacitance is unreal- 
istically large [2]. Golden and Halperin [7] obtain their 
results by mapping the two dot model to a single dot 
model. Here, we directly model the experimental set- 
up as a one-dimensional system and obtain the effective 



action exactly, by integrating out Gaussian degrees of 
freedom. Hence our results give the interaction energy 
directly in terms of the microscopic variables and ap- 
plied voltages. We find that in the weak tunneling limit, 
the peak splitting is proportional to the square-root of 
the tunneling conductance and in the strong tunneling 
limit, the deviation of the splitting factor from unity is 
proportional to the square-root of the deviation of the 
conductance (measured in units of 2e 2 /h) from unity. 
Model and effective action for coupled quantum dots 

Our model for a dot is essentially a short length of LL 
wire with barriers at either end, controlling the tunnel- 
ing to and from the dot. The two dot system is shown 
in Fig.(l), with three barriers in the LL wire between 
the two leads. The barriers are chosen to be (^-functions 
with strengths Gi and J. (Extended barriers do not 
change the results [7].) The electrons in the dot interact 
via a short range (Coulomb-like) interaction described 
by the Hamiltonian H = —vp J dx[^ L id x ipL — (L <-> 
R)] + g J dxp(x) 2 where p(x) = ip^ip is the electron den- 
sity and -0 = 4>Le~ lkFX + ^Re lkFX . Here, Vfl and ip L 
stands for fermion fields linearised about the left and 
right Fermi points. The electrons in the leads are free. 

The model is bosonised via the standard transfor- 
mation i\) v — rj^e 21 ^^" where v = R,L and 4> y are 
the bosonic fields. r\ v are the Klein factors that en- 
sure anti-commutation of the fermion fields. The La- 
grangian density for the bosonic field in 1+1 dimensions 
is £((/>; K,v) = {l/2Kv)(dt(f>) 2 -(v/2K)(d x (f>{x)) 2 , where 
K ~ (1 + g/irv F )- 1 ' 2 and (1 +g/irv F y/ 2 . The 

bosonised action is now given by [11-13] 

S = Jdr[S leads + S dots + S gates ], (1) 

with Si eads = (JZJ + J^)dxC((t);K L = l,v F ), 
Sdots = J° dl dxC(4>; K 1: v-i) + Jq 2 dxC((j); K 2 , v 2 ) + V + 
Jcos(20i : x) and S gat es = {gi/V^) f_ dl dx(x - 4>i) + 
{92/ V^) Jo* dx(4> 2 -x) , where V = S| =1 G l cos(2 + 
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(— l) l 2kFdi). Here Gi and J are the junction barriers 
defining the dots and fa = fa(—l) l di) and X = <j>(0) de- 
note the boson fields at those points. J controls tunneling 
between the dots and ranges from fully open (J = 0, sin- 
gle dot ) to fully closed (J = oo, two decoupled dots). 
The gate voltages gi control the density of electrons in 
each dot. 
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FIG. 1. Double dot system above shown schematically be- 
low with the dots modeled as LL with -ff;=i,2, between barri- 
ers Gi and J, and with two gate voltages gi. 

In the low T limit (T <C hvi/kBdi, hv 2 /kBd 2 ) all three 
barriers are seen coherently, and we expect the CB of the 
two dots to be coupled. The effective action of the system 
can be obtained by integrating out all degrees of freedom 
except those at the positions of the three junction barri- 
ers, following Ref. [11,12] and we find that 



Sef=S + J dT&U^iX- 

+V + Jcos(2V^x)} 



(2) 



where S a = £ Wn |w n |(# + ~4> 2 2 )/2{K L = 1) . The 
Fourier transformed tilde fields are defined by fa(r) = 
J2u e ~ luJ " T (t>i( UJ n)- The Ui = hvi/Kidi terms are the 
'mass' terms that suppress charge fluctuations on the 
dots and are responsible for the CB through the dots. 
Weak tunneling limit 

Here the barrier J is very large, and the field \ IS 
pinned at its minima - 



X = ny/ir = xd (n = integer). 



(3) 



A semi-classical expansion about these minima using x — 
Xci + Xq to second order in the fluctuations yields 

S ef ~ 5 + J dr^-gxfa + g 2 fa]/V^ + Xgidi - 92) /V* 

+-ZUU i [x 2 q + 2 Xq (Xci ~ fa) + (Xd - fa) 2 }/2 
+V-(2nJ) X 2 q } (4) 

where we have dropped all field- independent constants. 
We can now integrate out the quantum fluctuations in 
X , since we have only kept quadratic terms, and using 
J > U l ,g l , we get 

S ef ~ #0 + J dT[Y,j =1 l ^-{Xcl - faf + -^J^fo 

+{-gifa + g2fa)/V^ + V] . (5) 



Now, we define the following variables - rij = 
(-l) l+1 (Xci~fa)/V^+(kFd t )/ir and ft = ( Xc i + fa)/2 + 
(— l) I (fc_pcZi)/2 which can be interpreted as the charges 
and currents on the dots 1 and 2 respectively and in 
terms of which the action can be rewritten simply as 

Sef * So+ J dr^tiY^ - - Ei 2 nm 2 + V] (6) 
where 

E12 = U!U 2 /4J, k noi = {k F d i )/ir-g i /U i . (7) 

Ui = nUi plays the role of the charging energy or CB 
energy because in the absence of any mixing between r%\ 
and n 2l by tuning gi (equivalcntly n M ), the dot states 
with Hi and rij + 1 electrons can be made degenerate. 
This is the lifting of the CB for each individual dot. But 
the crucial term above is the term mixing m and n 2 . 
This tells us how the CB through one dot is affected by 
the charge on the other dot. The barrier terms can be 
more conveniently rewritten in terms of the total charge 
of the double dot and the current through the double dot 
system, when G\ = G 2 = G by defining the total charge 
and current fields as N — n\ + n 2 = (<p 2 — fa)/y/n + 
(k F L)/ii and 9 = (fa + fa)/2 + (k F l)/2^, where L = 
di + d 2 and I = d 2 — d\. With these redefinitions, the 
barrier term reduces to 



V = 2G cos^v 7 ?^) cos ttN 



(8) 



and So — > S' 



Y, Un K|(0 2 + (n/4)N 2 ). The action is 
now in a form where its symmetries are manifest. 

When the two dots are decoupled from each other 
(-E12 = or equivalently J — > 00), the action is sym- 
metric under the transformation 6^6 + ^/n, N — ► 
2iVo — N. This is because g\ and g 2 can be tuned so 
as to make n i, n 02 to be half- integers qi + 1/2, q 2 + 1/2 
and hence, No — noi + n§ 2 to be an integer. Thus the 
two CB's are lifted and the four charge states, (m, n 2 ) = 
{li, q 2 ), (qi + 1, q 2 ), (qi, q 2 + 1) (?i + 1, q 2 + 1) become de- 
generate and transport through the double dot system is 
unhindered. We can plot the points of maximum con- 
ductance through the double-dot system in the plane of 
the two gate voltages and as shown by the dotted lines in 
Fig. (2), we get a square grid. (The figure is for identical 
dots; in general, the grid is rectangular.) 

Now, we study the symmetries for a large finite J. 
In this case, the CB through each dot is influenced by 
interaction with the electrons in the other dot. Hence, 
we now find that the CB through the double dot system 
is lifted for two possible sets of gate voltages - 
• the states [(qi,q 2 ), (qi + 1,92)] and the states [(qi,q 2 ), 
(qi,Q2 + 1)] are degenerate simultaneously respectively 
when 



, 1 u 2 
noi=9l+ 2 + 8J <Z2 ' 



1 Ui 

n 02 = q 2 + - + —qi (9) 



This is one point where the CB is lifted and the conduc- 
tance is a maximum. 
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FIG. 2. Positions of the conductance maxima as a function 
of the gate voltages. The solid lines denote the splitting and 
the dot-dash line AB denotes the distance between two sets 
of split peaks. The splitting factor S is defined as twice the 
ratio of the former to the latter length. 

• Similarly, the states [(qi,q 2 + 1), (qi + 1,92 + 1)] and 
[(qi + l,q 2 ), (gi + 1, g2 + l)] are degenerate simultaneously 
when 

nai = qi + - + — (q 2 + 1), ^02 = 92+2 + gjWi + *) 

(10) 

This is another point where the CB is lifted and the con- 
ductance is a maximum. 

Note that the original resonance which occured when all 
four states were degenerate has now broken into two sep- 
arate resonances, at each of which only three states are 
degenerate. Plotting this in the plane of the gate volt- 
ages, we obtain a hexagon as shown by the dashed and 
solid lines in Fig. (2) in agreement with the experiments 
of Blick et al in [3]. (The figure shown is for identical 
dots; otherwise, the hexagons are not perfectly symmet- 
ric.) Thus, as a consequence of tunneling between the 
quantum dots, the conductance maxima splits into two. 

The splitting Sg of the maxima as a function of the 
tunneling barrier J (the solid line in Fig. (2)) can be com- 
puted from Eqs(9) and (10), and gives 

Sgi = 5g 2 = 71-LW2/8J => 8g = V2irUiU 2 /8J ■ (11) 

Hence, the splitting in gi and g 2 is identical and is in- 
versely proportional to the gate voltage controlling the 
tunneling barrier and directly proportional to the charg- 
ing energies of the two dots. However, when J becomes 
comparable to the other energy scales in the problem like 
Ui, the strong barrier analysis breaks down. 
Strong tunneling or weak barrier limit 

Here we assume that the barrier between the two dots 
is very small - i.e., d » Ui 3> J. So essentially, there is 
only one CB. However, we shall see that the spacing of 
the CB terms are affected by a non-zero J. We integrate 
out x ( m Eq.2) which is quadratic (when J is ignored) 
to obtain 

S eJ = EuJwnKfl 3 + \N 2 ) + U ef (N - N Q f/2 



+ 2G cos^y^) cos irN + J cosp^ + XirN + C] 
= S + V eJ . (12) 

Here, we have defined U ef = ttU 1 U 2 /{U 1 + U 2 ) = nU/2, 
X = (Ui — U 2 )HU X + U 2 ) = 0, N Q = k F L/ir - {g 1 jU 1 + 
92IU2) - k F L/Tr - g/U e f, C = 2{g 2 - .g 1 )/((7 1 + U 2 ) + 
(Ui + U 2 )k F l + (Ui - U 2 )k F L = 0. N and 6 have been 
defined earlier. C and N depend on the gate voltages 
and are the tunable parameters. The second equalities 
are for identical dots when U\ = U 2 = U and d\ = d 2 . 
gi = g 2 = g can be tuned by C. We shall use these values 
below for simplicity. 

Since, we continue to be in the strong barrier limit for 
the G term, the action has deep minima for integer TV 
(with appropriate 9). But the degeneracy between these 
minima is broken by the U e f and J terms leading to just 
one of them being preferred (as is the case for the usual 
CB without J). Just as earlier, the CB is lifted when 

V ef (6,N) = V ef (e + V^/2,N+l) . (13) 

When J = 0, this happens when Nq = half-integer. How- 
ever, when J 7^ 0, the above equality holds when 

7V = 1/2 - J/U ef , 3/2 + J/U ef , 5/2 - J/U ef (14) 

which implies that the splitting 5g is given by 

Sg = U ef - 2J . (15) 

It is easy to see from Eq.(14) that periodicity is restored 
when two electrons are added to the system, but peri- 
odicity when a single electron is added to the system is 
broken by the J term. This tells us that the main ef- 
fect of a small barrier within a single dot, is to change 
the spacing of the peaks with two peaks coming closer to 
each other and pairs of peaks receding from each other. 
But the distance between pairs of peaks remains 2U. As 
J increases towards 00, the peaks which come closer to 
one another slowly merge into one, and we are left with 
half as many peaks, which is what one would expect when 
the dot size is halved. 
Discussions and conclusions 

Thus, both from the weak and strong tunneling lim- 
its, we get a consistent picture. The splitting is zero for 
J = 00 and maximum when the tunneling goes to one 
(J = 0) and is just the spacing of a single dot of length 
d\+d 2 . As a function of the tunneling barrier, the split- 
ting saturates linearly (proportional to J) for low barri- 
ers (strong tunneling) and goes as 1/ J for weak tunneling 
(strong barriers). (Our model does not include inter-dot 
capacitance since it is expected to be very small for the 
experiment [2].) In terms of the barrier conductance G&, 
simple quantum mechanics shows that (in units of 
2e 2 /h) falls off from unity as J 2 in the low barrier limit 
and increases from zero as 1/ J 2 in the weak tunneling 
limit. So in the weak tunneling limit, using Eq.(ll) and 
Fig. 2, the splitting factor S = Sg/(U + Sg) - Sg/U (for 
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Sg <C U) oc 1/J oc VG&- In the strong tunneling limit, 
S = 2Sg/2U ef = 1 - 2J/U ef . So 1 - S oc J oc VI - G b 
Thus, in both limits, we relate the splitting factor to the 
barrier conductance. This is the central result of this 
letter. 

In Fig. 3, we compare our theoretical prediction of S 
with the experimental curve. In the weak tunneling limit, 
a one-parameter \ 2 fit of our prediction to the lowest five 
points (upto S ~ 0.3) gives a goodness of fit of 84% . In 
contrast, a linear fit to the same data ( predicted in Ref. 
[7]), gives a goodness of fit of only 75%. Similarly, in 
the strong tunneling limit (S above 0.7), a one parame- 
ter fit of our prediction (which, in fact, is in agreement 
with the scaling analysis of Ref. [14]) gives a goodness 
of fit of 96%, again, considerably better than the 81% 
for the logarithmic prediction in Ref. [7]. However, al- 
though our modelling is more realistic than earlier ones 
and the agreement of our predictions with the experimen- 
tal data, both in the strong and weak tunneling limits is 
impressive, better data is required for a conclusive proof. 
Further, we predict that for the weak tunneling case, the 
interdot interaction energy is directly proportional to the 
charging energies of the dots and hence inversely propor- 
tional to the sizes or capacitances of each of the individual 
dots. 

The charging energies and the Fermi velocity obtained 
from Ref. [2] are ~ AOO^eV. and (v F ) ld ~ 2.3xl0 7 cm. 
Thus using the relation, U = hv/Kd, for a single dot, 
we find the LL parameters, Vi = v% ~ 3.1 X 10 7 cm/ sec 
and Ki — K 2 ~ .74 for identical dots. In the weak 
barrier limit, where the splitting is small, we find that 
E12 (in Eq.7) ranges from ~ 4 — AO^ieV, which is roughly 
(10~ 2 — lO^ 1 )?/. The LL parameter values remain almost 
unchanged in the range of gate voltages used in the ex- 
periments. The values can be confirmed by studying the 
conductance through the double dot system in the 'high' 
temperature limit, T Td — hv/ksd ~ 0.6K , when the 
electron transport through each of the dots is incoher- 
ent. The conductance should scale as T 2i - X I K ~^ ~ T o m 
in the weak tunneling case and as T 2 ( K ~^ ~ j"-o.5i m 
the weak barrier case. LL behaviour can also be probed 
if the conductances are studied as a function of the size 
of the dots at 'low' temperatures T <^Td since LL the- 
ory predicts explicit size dependent power laws at low 
temperatures. 

In conclusion, in this letter, we have obtained the effec- 
tive action of a system of coupled quantum dots, in both 
the weak and strong tunneling limits. We have shown 
that in the presence of inter-dot tunneling, the peaks de- 
noting the conductance maxima split. The split is maxi- 
mum when the tunneling between the dots is maximum 
- i.e., when there is no barrier at all between the dots. 
We have also computed the splitting as a function of the 
gate voltage controlling the tunneling between the dots 
in both the weak and strong tunneling limits, and find 
that the splitting factor is proportional to the square-root 
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FIG. 3. Splitting factor S as function of the gate voltage 
controlling tunneling. The data points are from Ref.(2). The 
solid line denotes the square-root of the measured conduc- 
tance (multiplied by scaling factor 0.56 for best fit), and the 
dotted line just denotes the measured conductance (multiplied 
by 0.60 for best fit). 

of the conductance in the weak tunneling limit and the 
deviation of the splitting from unity is proportional to 
the square-root of the deviation of the conductance from 
unity in the strong tunneling limit. This agrees with the 
experimental data on double dot systems. We have also 
extracted the Luttingcr parameters from the experiment 
and have predicted temperature power laws for the same 
set-up which can be experimentally tested. 
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